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Abstract 

We demonstrate numerically that for Ising spins on square lattices with 
ferromagnetic second neighbour interactions and random near neighbour in- 
teractions, two dimensional Ising spin glass order with a non-zero freezing 
temperature can occur. We compare some of the physical properties of these 
spin glasses with those of standard spin glasses in higher dimensions. 
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In dimension two, an Ising spin glass (ISG) with random near neighbour interactions 
does not order above zero temperature Implicitly this has been taken to mean that 

two dimensional ISGs with ordering at finite temperatures do not exist. Here we show that 
this assumption is unfounded by exhibiting a family of 2d Ising systems with a different 
set of interactions, which show spin glass ordering at finite temperatures. We find that 
the physical properties of these 2d ISG systems ressemble closely those of standard ISGs in 
higher dimensions. 

Consider a square lattice of Ising spins with second nearest neighbour ferromagnetic 
interactions of strength J. This clearly consists of two inter-penetrating but non-interacting 
square sublattices which will each order ferromagnetically at the Onsager Curie temperature 
T c = 2.27 J. Now introduce random near neighbour interactions of strength ±AJ. At low 
temperatures when the spins are frozen, each sublattice will exert effective random fields on 
the spins of the other sublattice. It is well established |3[ that a 2d Ising ferromagnet breaks 
up into finite size domains under the influence of a random field, however small. At low 
temperatures each sublattice of the present system will then be ordered in randomly arranged 
ferromagnetic domains of finite size. However in contrast to the true 2d random field Ising 
system, here the "effective random fields" will average to zero at high temperatures, so we can 
expect a critical temperature T g separating a high temperature paramagnetic regime (with 
the autocorrelation function q(t) relaxing to zero at long time t) from a low temperature 
regime where q(t) does not go to zero, corresponding to ordering with a non-zero Edwards- 
Anderson order parameter. Below we will show numerical evidence that for small and 
moderate A, T g is non-zero. In the low temperature state there is no long range ferromagnetic 
order so the system is certainly not a ferromagnet; we will choose to call the frozen state a 
"spin glass" , although random systems with short range ferromagnetic or antiferromagnetic 
order are often referred to as "cluster glasses" in an experimental context. 

One of the established numerical techniques for determining the ordering temperature 
in ferromagnets and spin glasses is through finite size scaling using the Binder cumulant [|J . 
For sets of samples of different sizes L, the fluctuations of the autocorrelation function q(t) 
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in equilibrium are recorded. The Binder cumulant at temperature T is denned as: 



9l(T) 




) 



(1) 



and is dimensionless; for sets of samples with fixed L values, the series of curves gi{T) 
as a function of temperature intersect at the ordering temperature. We have carried out 
simulations using heat bath dynamics and sequential spin by spin updates on samples up 
to L — 12. 500 samples were used at each size for each value of A and standard precautions 
were taken to assure that thermal equilibrium was achieved @]. We show in figure [l] the 
gL data for three values of A and it can be seen that for the smallest value there is a clear 
intersection point. For A = 0.7 the gL curves come together without any clear fanning 
out at lower temperatures, but we consider that the meeting point represents the freezing 
temperature here too. For A = 1 it is hard to say if the curves will intersect above T = 0. 
(Very similar behaviour with almost unobservable fanning out below T g is seen in 3d ISGs 
The behaviour of the freezing temperature as a function of A is shown in the inset of 
figure |2|; at small A, T g extrapolates back to the Onsager T c , and with increasing A, T g drops 
as might be expected, because for high enough A we will recover the 2d random interaction 
near neighbour ISG with its zero temperature ordering. It appears that the limiting value 
of A for non-zero T g is in the region of 1.0. 

In conventional second order transitions gL{T c ) has a universal value for all systems in 
a given class; for instance for the 2d Ising ferromagnet the magnetism Binder parameter 
9Lm(T c ) = 0.91603 || and we find the autocorrelation Binder parameter (7l(T c ) = 0.42. We 
can note in figure ([[]) that giiTg) changes with A so for the present systems of 2d ISGs 
universality is not obeyed. For the range of A that we have studied, ghiTg) appears to be 
tending regularly to a value near 1 as T g is tending to zero. At still higher A values ^l(0) 
will presumably tend to 0.83, which is the standard 2d ±J ISG zero temperature value |§f§]. 
We can note that a closely related non-random 2d Ising system, having first and second 
neighbour interactions with regular frustation, shows non-universality of critical exponents 
when a control parameter representing the ratio of the two interaction strengths is varied 
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||. There is also strong evidence for non-universality of the critical exponents in canonical 
ISGs at dimension 3 and above JF|. We have not yet estimated the values of the critical 
exponents for the present family of 2d ISGs. 

Once the existence of finite ordering temperatures established, we can carry out simu- 
lations to measure the standard physical properties for these 2d ISGs in order to compare 
with those of other ISGs. In figure (§) we show the specific heats for the same three values 
of A. It can be seen that the specific heats show broad maxima at temperatures higher 
than the respective ordering temperatures, and smooth featureless behaviour through the 
freezing transition. This is precisely what is observed for standard ISGs in dimension 3 |§, 4 
and 5 || and appears to be a universal characteristic of the ISG transition below the upper 
critical dimension. For the present systems the maximum of the specific heat corresponds 
to a saturation of the ferromagnetic correlation length in each sublattice as the sample is 
cooled; the ordering process can be thought of as consisting of first the formation of finite 
size ferromagnetic domains, and then of a freezing at a lower temperature where the domain 
movements are hindered by the random near neighbour interactions. 

Another parameter accessible from simulations is the damage spreading |1U|| . Damage 



spreading is measured numerically by using heat bath dynamics and applying the same 
random update parameter at each update step to two replicas of the same system; D(T) 
is defined as the normalised Hamming distance between the two replicas (i.e. the fraction 
of spins having opposite orientations) at long times when this protocol is applied. While 
ferromagnets have D(T) = down to T c , standard ISGs have been found to show non-zero 
D(T) until temperatures well above T g | IoyiI |; it has been pointed out that for ISGs in 



dimension 3 and 4, as the temperature is reduced D{T) tends to 1/2 at a temperature very 
close to T g [jO]]. We have studied damage spreading in the present 2d ISGs. Two random 
initial replicas of each 2d system were chosen; these replicas were annealed independently 
to thermal equilibrium at each temperature, and then the damage spreading procedure was 
applied for long times. The results, figure (^|), show that the damage spreading temperature 
Tjj where D(T) first becomes non-zero is much higher than T g , and that in each case on 



cooling D{T) tends to 1/2. For A = 0.5 this temperature is indistinguishable from the Tg 
we have estimated from the Binder cumulant data. D{T) then remains equal to 1/2 for 
lower temperatures. The behaviour is similar but less clear cut for the other two values of 
A (the lowest temperature points may be slightly low because of the very long times needed 
to achieve complete thermal equilibrium). The general behaviour of the damage spreading 
parameter appears to be the same in these 2d ISGs as in ISGs in dimensions 3 and 4. The 
use of the empirical criterion D(T g ) = 1/2 gives a convenient method to check T g values. 

An advantage of the present systems over conventional ISGs is that "snapshots" of the 
instantaneous configurations of the spins are informative because we are in dimension 2. We 
show in figure (f|) two snapshots of the same A = 0.5 system at temperature T = 1.5, well 
below T g . The images were produced by slowly cooling two independent random replicas of 
the same system by successive steps to the final temperature. As could be expected there 
are four types of domains : if we call one sublattice A and the other B then the domains are 
A up B up (white), A up B down (white/black checkerboard), A down B up (black/white) 
and A down B down (black). Here the domain sizes are typically 25 by 25 spins. (The 
simulations confirm that the typical domain size below Tg gets smaller as A is increased). 
For this sample at this temperature the two configurations are almost frozen; averaging 
over a further 10000 Monte Carlo Steps per spin (MCS) only blurs the domain frontiers 
somewhat. To a good approximation the images then represent two Gibbs states of the 
system. By inspection it can be seen that the two states are neither quasi-identical, nor are 
they quasi-mirror images of each other with all the spins reversed. The detailed structures 
of the two states are quite different suggesting that there are an infinity (or at least a large 
number) of alternative Gibbs states for this ISG. This conclusion appears to be incompatible 
with the two ground state image which has been proposed heuristically for ISGs E2[ but the 



phase space could ressemble that of the Parisi solution of the infinite range SK model [To 

In conclusion, we have demonstrated that 2d ISG systems with non-zero freezing tem- 
peratures exist. We have shown data for a certain number of parameters for a new family 
of 2d ISGs where we have found that the observed behaviour is very similar to that seen 



in simulations on standard ISGs in dimension 3. A more comprehensive study of the 2d 
ISGs and a detailed comparison with results from higher dimensions should establish what 
properties are always necessarily associated with spin glass ordering and should help under- 
standing of the spin glass phenomenon in general. As far as experimental comparisons are 
concerned, it might be worth reconsidering the data on 2d ISG materials without making 
the a priori assumption that the ordering temperature is necessarily zero. 

Finally, the general philosophy behind the present study may help to resolve a long 
standing paradox. Real three dimensional spin glass materials where the spins can be con- 
sidered to be basically Heisenberg show well defined non-zero freezing temperatures; on the 
other hand in numerical work and on theoretical grounds it appears that Heisenberg spins 
with random near neighbour interactions should not show a freezing transition at a finite 
temperature |L4j (but it should be noted that with vector spins chirality should also play a 



role |15j and that a simulation of a Heisenberg spin glass with RKKY interactions shows a 
finite T g |16[]). Now we have just seen that if we have a set of random near neighbour inter- 
actions only, a 2d Ising system has no finite temperature ordering, but the system can be a 
spin glass with a finite T g if there are both random and non-random interactions. Perhaps 
the same rule applies to 3d Heisenberg systems. It may well be that all real life Heisenberg 
spin glasses have both random and non-random interactions (and so in a sense might be 



classified as "cluster glasses" ||17|| ). Thus even in the canonical CuMn alloys, the short range 



magnetic ordering which has been observed with neutrons [18j may be essential for the spin 
glass transition to occur at a finite temperature. 
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FIGURES 

FIG. 1. The measured Binder Cumulants as a function of size L and temperature T for 
A = 0.5,0.7 and 1.0 respectively. 

FIG. 2. The specific heat averaged over samples with L = 40 for A = 0.5, 0.7 and 1.0 (circles, 
squares and triangles respectively). 
Inset T g as a function of A 

FIG. 3. The damage spreading parameter D(T) as a function of temperature (see text) for A 
=0.5 (triangles), 0.7 (circles) and 1.0 (squares) for 40 x 40 spins. 

FIG. 4. Snapshots of two replicas A,B of the same 100 x 100 spin system with A = 0.5 at T=1.5 
after 10000 mcs anneal. The two replicas were cooled independently from two initially randomly 
chosen configurations On each site black is up white is down. The third frame is the difference 
between the two configurations A and B; black: Sf = S- 3 , white Sf = —S- 3 . 



9 



1 .0 I ■ 1 ■ 1 ■ r 




0.0 1 ' 1 ' ' ' ' ' 1 

1.0 1.5 2.0 2.5 3.0 

Temperature 



1 .0 I ■ 1 ■ 1 ■ r 




Temperature 



1 .0 I ' 1 1 1 ■ 1 1 r 




Temperature 



